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Abstract. We present a model independent study of the chiral condensate evolution in a hadronic gas, in 
terms of temperature and baryon chemical potential. The meson-meson interactions are described within 
Chiral Perturbation Theory and the pion-nucleon interaction by means of Heavy Baryon Chiral Pertur- 
bation Theory, both at one loop. Together with the virial expansion, this provides a model independent 
systematic expansion at low temperatures and chemical potentials, which includes the physical quark 
masses. This can serve as a guideline for further studies on the lattice. We also obtain estimates of the 
critical line of temperature and chemical potential where the chiral condensate melts, which systematically 
lie somewhat higher than recent lattice calculations but are consistent with several hadronic models. 
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On a recent work [11 we approach the question of the 
phase diagram of QCD. In particular, we study the tran- 
sition from the hadronic phase, in which the chiral sym- 
metry is spontaneously broken, to a phase in which it is 
restored. This is done within Chiral Perturbation Theory 
(ChPT), the low energy effective theory of QCD. We also 
gave estimates for the melting curve of the quark conden- 
sate in presence of a baryon chemical potential. 

Let us briefly describe ChPT. The spontaneous chi- 
ral symmetry breaking of QCD requires the existence of 
eight massless Goldstone Bosons (GB), that can be iden- 
tified with the pions, kaons and eta. These are thus the 
most relevant degrees of freedom at low energies. In ad- 
dition, there is an explicit symmetry breaking due to the 
non- vanishing quark masses that give rise to a small mass 
for the pions, kaons and eta, which are thus just pseudo- 
GB. ChPT is built as the most general derivative and 
mass expansion over the spontaneous symmetry breaking 
(SSB) scale = dTr/^r ~ 1.2 GeV which respects the 
symmetry constraints of QCD. Baryons can be included 
in this formalism, although its treatment is more involved 
due to their large masses. Within Heavy Baryon ChPT 
(HBChPT) |2| this problem is overcome for the meson- 
baryon interaction by an additional expansion on inverse 
powers of the baryon mass. This approach has proven very 
successful, and works remarkably well within the meson 
sector, and with a somewhat slower convergence in the 
meson-baryon sector. 

Next, the thermodynamics of the hadron gas is de- 
scribed by the virial expansion 4,5J. This is a low density 
series expansion for the grand canonical potential of a rel- 
ativistic interacting multicomponent gas. We will assume 
that only the strong interactions are relevant, and that the 



baryon density defined as ub — is conserved. The ther- 
modynamics of the gas will thus depend on the tempera- 
ture T (usually written in terms of its inverse P ~ I /T) 
and a baryon chemical potential /i^. The virial expansion, 
written in terms of the pressure, reads: 
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where — exp/3{fj.i — Mi), with Mi the mass of the i- 
th species, /i^ — 0, ±/is for mesons/ (anti)baryons respec- 
tively. The free virial coefficients 
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, with rji — ±1 for bosons/mesons and gi the degeneracy of 
the i-th species, correspond simply to the series expansion 
of the pressure for a free gas, while the interactions are 
taken into account in the term: 



^0(Mi+Mj) poo 



dEE'^Ki{PE)A'^{E), (3) 



Mi+Mj 



where ^^^1(0;) is the modified Bessel function (of the sec- 
ond kind), and A'^ = T.i,.j,si'^I + 1)(2^ + 1)^*^' are the 
ij — > ij elastic scattering phase shifts for a state ij with 
well defined isospin /, total angular momentum J and 
strangeness 5, defined so that (5 = at threshold. 

The dependence of the chiral condensate with temper- 
ature and chemical potential is written in terms of the 
pressure of the gas as follows: 



{qq) = (0|gg|0) 1 + E 
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Fig. 1. Crude estimate of the region of validity of the virial 
expansion for a free gas. We have excluded the regions in which 
the virial expansion breaks, but also those regions in which 
nuclear density is large enough so that higher order effects are 
relevant. 



Fig. 2. Melting line in the {fj,B,T) plane of the chiral conden- 
sate for an SU(2) gas of hadrons. We show the result for the 
free gas, but also adding the tttt ChPT SU(2) interaction to one 
loop, and the line resulting from adding the nN interaction to 
third order in HBChPT. 



where the constant F, which is the pion decay constant 
in the chiral limit, is introduced for further convenience, 

and the coefficients Ch = —F'^^^^{0\qq\0)~^ encode the 
hadron mass dependence on the quark mass. Numerically, 
these coefficients amount to c^r — 0.9lg 4, ck = O.StQ;^, 
Cri = 0.4^g'7, for mesons [6 , and = 3.6t.\'g for nuclc- 
ons [1]. 

First of all, we want to estimate the region of valid- 
ity of our calculations. For this, we check the points for 
which the virial expansion breaks down, by computing the 
difference between the free condensate calculated exactly 
and with the virial expansion. In Figure [1] we show, as a 
dashed and a dashed-dotted lines, the points where the 
virial expansion is off the exact calculation by 5% and 
20%, respectively. The expansion rapidly deteriorates af- 
ter that. Also, as we want to neglect NN interactions, 
we have to restrict ourselves to the region where nucleon 
density is small enough. In ref. [7] it is shown that, for 
densities below po/2, with po — 0.16 fm^' , and as far as 
we are only concerned with the quark condensate, the nN 
interaction dominates over the NN interaction. Thus, we 
restrict the validity of our calculation to the points be- 
low the po/2 line. This line is shown in the Figure as a 
dotted line, together with the po line, which is shown for 
illustration as a continuous line. The white area in Fig. 2 
corresponds to our estimated "validity region" . 

In Figure [2] we show in the {(Xb^T) plane the con- 
densate melting line of the free SU{2) gas composed of 
pions and nucleons as a dash-dot-dotted line. We have 
then added the effects of the tttt interaction, which yields 
the dashed line of the Figure. Finally, the ttN interaction 
was included. We see that the effect of the nN interac- 
tion is small, but noticeable (~ 6 MeV at /is = 0) at 
low baryon chemical potentials, and maximum at around 



Pb — 600 — 700 MeV, where it produces a decrease in the 
melting temperature of about 40 MeV with respect to the 
gas without -kN interactions. Up to a chemical potential 
oi fj,B = 40 — 50 MeV, the region of relevance for Rela- 
tivistic Heavy Ion Collisions, the decrease in the melting 
temperature amounts roughly to 10 MeV when we include 
the irN interaction. It is important to note that, although 
the critical lines are actually extrapolations, at low tem- 
perature and chemical potential our calculation is model 
independent, and should be quite accurate. Still, we show 
the melting lines since they are useful to visualize and help 
quantifying the relative size of each contribution we add 
into the gas. 

In a real gas, we should consider all hadrons. We will 
do this by including them as free particles. The only ex- 
ception to this will be the kaons and etas, which are abun- 
dant enough up to 200 MeV to deserve a separate treat- 
ment and include their interaction with a pion. All other 
interactions are suppressed by Boltzmann (thermal) and 
Ch/Mfi factors, and are therefore not included in our treat- 
ment. At very large fj,B the heavier nucleons may not be 
Boltzmann suppressed, and their interactions may become 
important, but this is outside the scope of this work. In 
Figure[3]we show the two lines corresponding to the SU{2) 
gas with TTTT interactions (dashed) and with tttt and ttN in- 
teractions (continuous). We then plot the line correspond- 
ing to a gas in which also kaons, etas, and their interac- 
tions with a pion, now in SU (3), are included (dotted line). 
Note that the effect is pronounced at ps = 0, with a de- 
crease of the melting temperature of around ~ 13 MeV, 
but becomes negligible at very high baryon chemical po- 
tential. Finally, we have added the effects of other heav- 
ier hadrons, included as free particles, since they suffer 
a larger Boltzmann suppression. We have estimated that, 
for heavier hadrons, dMh/drh ~ ctN^ ^, where a is an adi- 
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Fig. 3. Our final result for the chiral condensate melting 
line in the {fiB,T) plane. Starting from the SU(2) gas with 
TTTT interactions, we show the effect of nN interactions and of 
adding kaons and etas. The dark area covers the uncertainty 
due to heavier hadrons estimated as explained in the text. 



mensional constant and ^ is the number of valence u or 
d quarks in the hadron. Based on the values of a that we 
have calculated for the pions, kaons, etas and nucleons, we 
have estimated that, for heavier hadrons, a ~ 0.5 — 2.5. 
The band in Figure [3] corresponds to this range of val- 
ues for a. In order to properly analyse the uncertainties 
in our calculation, we also have to include the uncertain- 
ties that come from our imperfect knowledge of the chiral 
parameters. This was done by adding in quadrature the 
errors coming from independently varying each of the pa- 
rameters. The final results are plotted in Figure 21 where 
the dark band corresponds to the uncertainties coming 
from the a parameter in heavy hadrons and the dashed 
line includes also the uncertainty coming from the chi- 
ral parameters. Note that the error is highly asymmetric: 
even though every contribution lowers the melting tem- 
perature, the melting of the condensate accelerates near 
the melting point, and thus even if the contribution has a 
similar size to the previous one, its effect on the conden- 
sate melting seems smaller. 

It is interesting to compare our results with those ob- 
tained by using lattice calculations. In Figure [5] we show 
our curves with their errors, plotted against other results 
recently appeared in the literature. Our approach agrees, 
within errors, with the hadron resonance gas model de- 
scribed in [5]. However, ChPT calculations [MSJi includ- 
ing ours [T], yield melting temperatures that are system- 
atically above the results from the lattice calculations [9] . 
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Fig. 4. Analysis of uncertainties. To the dark area, which 
corresponds to the case described in Fig. O we have added 
the uncertainties coming from the chiral parameters, shown as 
dashed lines (see ^1 for discussion). 
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Fig. 5. Comparison between our results and others in the lit- 
erature. Our melting temperatures are consistent within errors 
with ref. [S] but systematically higher than [9]. 
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